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1 Introduction 



Nuclear spins may exhibit several coherent phenomena, occurring at nuclear mag- 
netic resonance frequencies, which have their direct counterpats in resonant atomic 
systems"*^ . For instance, nuclear free induction is an analog of atomic free induction 
and the spin echo is an analog of the photon echo. These analogies are due to the 
fact that an ensemble of identical spins forms a collective of finite-level objects, sim- 
ilar to a resonant system of atoms or molecules. Such a finite-level nonequilibrium 
system can, under special conditions, behave as a collection of coherent radiators. 

Similarly to the case of resonant atoms, one may distinguish two principally dif- 
ferent ways of dealing with spin assemblies. One situation would be when spins are 
near their equilibrium state, with a nonequilibrium perturbation produced by a res- 
onant alternating field. This is a typical situation of nuclear magnetic resonance . 
Another possibility could be if spins are initially prepared in a strongly nonequilib- 
rium state, e.g. being polarized against a constant external magnetic field. In the 
latter case, the polarized spins resemble a system of inverted resonant atoms. 

One of the most interesting effects exhibited by a system of inverted atoms is 
superradiance, when the radiation intensity is approximately proportional to the 
number of atoms squared, while the duration of a superradiant pulse is inversely 
proportional to this number. The possibility of atomic superradiance was predicted 
by Dicke , and nowadays it is well studied both theoretically and experimentally. 
There exists a vast literature on the subject, whose description can be found in the 
recent books^'^. A natural question that arises is: If atomic and spin systems are 
so similar, then could a kind of spin superradiance be realized? 

For spins to behave coherently, there must exist a cause correlating their motion. 
In the case of atoms, such a correlating mechanism is caused by the photon exchange 
through the common radiation field. This exchange results in the formation of ef- 
fective interatomic correlations collectivizing atomic radiation. For spins, however, 
their magneto-dipole radiation is too weak to develop noticeable spin correlations. 
This concerns both nuclear as well as electron spins. Such photon-exchange interac- 
tions are negligible as compared to disordering dipolar spin interactions. How then 
might the spin motion be collectivized? 

It is, of course, possible to force an ensemble of spins to develop coherence by 
imposing an initial condition of longitudinal magnetization, and using an rf pulse 
to produce a macroscopic transverse magnetization. But this would result in the 
standard free nuclear induction, with the coherence being lost during the dephasing 
time T2. However, free nuclear induction, although being a coherent process, is not 
superradiance with one of the main characteristic features of a short radiation time 
of a superradiant pulse Tp <^ T2, shorter than the transverse dephasing time. To 
make the superradiant pulse time Tp so short, some internal nonlinear correlating 
mechanism has to be involved. 

Mutual spin correlations can arise owing to a feedback field formed by a reso- 
nant electric circuit coupled to the spin system and tuned to the Zeeman transition 
frequency of spins. The role of such a coupling in magnetic resonance experiments 
has been analysed by Bloembergen and Pound^. They showed that, in the presence 
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of an electric resonant circuit, coupled to a spin system, the signal of nuclear in- 
duction can be damped in a time much smaller than T2. Since this shortening of 
the induction damping time is due to collective effects, caused by the spin coupling 
through the resonator feedback field, the resulting process may be termed collective 
induction. This effect has been observed in a number of substances, as has been 
reviewed . The process already possesses one of the prerequisits of superradiance, 
i.e. a short radiation time Tp < T2. Usually, the induction signal starts at t = 0, hav- 
ing there its maximum intensity, while a superradiant pulse is always separated from 
the time origin by a delay time to > 0. The signal of collective induction can also 
be peaked at a delay time well separated from t — 0. However, there is a principal 
difference between collective induction and superradiance: In collective induction, 
coherence is induced by external sources, while in superradiance coherence develops 
in a self-organized way owing to internal correlations. The self- organized, sponta- 
neous, nature of superradiance is one of its most important features ^ . In general, 
one includes into the class of superradiant the processes that are triggered by initial 
external pulses, but with a compulsory requirement that the imposed initial coher- 
ence be very weak, playing just the role of a trigger. Thus, collective induction, 
though being a collective coherent process, is not yet superradiance. 

The principal characteristics of superradiance, which apply both to atomic as 
well as to spin systems, can be summarized as follows: 

1. Superradiance is collective, coherent radiation by an ensemble of radiators. 

2. It is a spontaneous process developing in a self-organized way. 

3. The maximal intensity of a superradiant burst is proportional to the number 
of radiators in the power larger than one. 

4. The duration of a superradiant pulse is essentially shorter than the dephasing 
tine T2. 

5. A superradiant pulse has a peak at a finite delay time. 

In atomic systems, one may distinguish different types of superradiance, all hav- 
ing the same characteristic features, including a short pulse time Tp < T2 and finite 
delay time to > 0. First of all, superradiance can occur as a transient process or 
as a lasting repeated effect, depending on whether there is no external nonresonant 
pumping permanently applied to the system or, respectively, there exists such a 
pumping supporting atomic inversion. In the first case, the system is prepared in 
an inverted state, after which no nonresonant pumping is involved. Then, if the 
relaxation process is gently promoted at t = by an external pulse imposing weak 
initial coherence on the system, triggered superradiance may develop. When no ini- 
tial coherence is thrust on the system, but superradiance appears as a completely 
self-organized effect starting from a purely incoherent state, then this is named pure 
superradiance. In the case, when the system is subject to the action of a nonresonant 
pumping constantly supporting atomic inversion, the regime of pulsing superradiance 
may arise, with a long series of superradiant bursts. 
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It has been of great interest to find out if a kind of spin superradiance, occurring 
in spin systems similarly to this phenomenon in atomic assemblies, could be realized 
experimentally and, if so, how should it be described theoretically. Despite several 
similarities between spin and atomic ensembles, they are, nevertheless, rather dif- 
ferent. For example, spontaneous radiation, starting the process of pure atomic 
superradiance, is absent in spin assemblies. Therefore, one of the most intriguing 
questions has been if pure spin superradiance can exist and, if so, what would be its 
origin. 

The phenomenon of spin superradiance, being realized, could be of interest by its 
own and, in addition, it could be employed for several applications, as is discussed 
in section 11. For example, these applications could include: 

(i) Investigation of materials characteristics by measuring the relaxation param- 
eters that are specific for spin superradiance and do not arise in other types of spin 
relaxation. 

(ii) Fast repolarization of targets used in various scattering experiments of high- 
energy physics. 

(iii) Construction of spin masers producing coherent radiation at radiofrequen- 
cies. 

(iv) Creation of sensitive detectors of weak external pulses, based on the mech- 
anism of triggered spin superradiance. 

(v) Method of information processing, derived from the feasibility of regulat- 
ing the number of and the intervals between superradiant bursts in the regime of 
punctuated spin superradiance. 

2 Experimental Observation 

Q 1 A 

Historically, the regime of pulsing spin superradiance was observed first°"'". This 
was done for a ruby crystal AI2O3, with the Cr^"*" paramagnetic admixture. The 
active nuclei were 27A1, with spins / = 5/2. If the ruby crystal is oriented in an 
external magnetic field such that a fully resolved quadrupole structure of its five 
Am = ±1 NMR transitions can be observed, and a resonant circuit is tuned to a 
selected NMR line, then ^^Al spins form a fictitious two- level system. In experiments, 
the resonant circuit was tuned to the central line. The density of ^^Al 

nuclei was pAi = 4.4 x 10^^ cm~^ and that of Cr^+, pcr = 8.6 x 10^^ cm~^. The 
measurements were performed in a static magnetic field of about Bq = 1.1 T, which 
for the ^'''Al magnetogyric ratio 7 x 10''' s~-^T~^ makes the NMR frequency 8 x 10''' 
Hz. The temperature range was from 1.6 to 4.2 K. The ringing time of the resonant 
circuit was r = 10~^ s, and the quality factor Q was varied between 60 and 200. 
The coil filling factor was 77 = 0.55. The transverse spin damping time was T2 = 

3 X 10~^ s. The spin population inversion of the ^''Al nuclear spins was achieved by 
dynamic nuclear polarization, by means of powerful microwave radiation supplied 
to the sample in the vicinity of a selected Cr^+ electron spin resonance line. The 
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corresponding longitudinal spin pumping time was Tj* = 0.1 s. The emission of a 
long train of superradiant bursts was observed, with the delay time of the first pulse 
being = 2 x 10^'^ s and its duration being = 2.6 x 10~^ s in the case of the 
quahty factor Q = 60, while = 8 x 10"^ s for g = 200. 

11 12 

The first observation of pure spin superradiance was accomplished in Dubna ' . 
A dielectric propanediol C3H8O2, with a paramagnetic admixture of Cr^+, was used 
as an active substance. This material possesses a high concentration of protons, with 
density pn = 4: x 10^^ cm^^. The proton spins J = 1/2 played the role of radiators. 
The admixture of Cr^+, with the density pcr = 1-8 x lO^'^ cm~^, was employed 
for the purpose of dynamic nuclear polarization of proton spins. Experiments with 
the same material propanediol were repeated in Saint Petersburg"*^ . In all these 
experiments , electric circuits were used as resonators, with quality factors Q 
between 100 and 600. The filling factor was rj = 0.6. The volume of samples varied 
from 0.5 cm^ to 12 cm^. For external static magnetic fields Bq in the interval between 
0.5 T and 2.64 T, with the proton magnetogyric ratio 2.675 x 10^ s-^T-^ the NMR 
proton frequencies were between 1.3 x 10^ Hz and 7.1 x 10^ Hz. The experiments were 
carried out at low temperatures of 0.05 K to 0.1 K. The inverted proton polarization 
reached 90%. Cooling of the sample resulted in strong suppression of the nuclear 
spin-lattice relaxation. The related longitudinal relaxation time was Ti = 1.8 x 10^ 
s at So = 0.5 T and Ti = 1.8 x 10^ s at So = 2.64 T. The transverse dephasing 
time was T2 = 0.85 x 10~^ s. After preparing an inverted sample, with the proton 
polarization directed against an external static magnetic field Bq, the latter was 
scanned with the velocity about 5 x 10~^ Ts~^. When in the process of this slow 
adiabatic scanning the field value Bq reached that one for which the proton NMR 
frequency coincided with the circuit natural frequency, the resonance condition was 
met. Then a powerful superradiant burst was produced, with the duration time 
about Tp = 1.3 X 10"*" s. 

Similar experiments, observing pure spin superradiance, were accomplished in 
Bonn^^, with the target materials butanol C4H9OH and ammonia NH3. These 
materials are also rich with protons, with density pj^ fa 3 x 10^^ cm~^. By means 
of dynamic nuclear polarization, proton polarizations of up to 99% were achieved. 
The samples were cooled to low temperatures, the spin-lattice relaxation time was 
Ti ^ 3.6 X 10^ s for protons in butanol and Ti ^ 10^ s for protons in ammonia. The 
resonant electric circuit, with the quality factor Q = 33, had a resonance frequency 
of 1.6 X 10^ Hz. 

In different experiments, slightly different setups were employed. A typical 
arrangement ^^'^^ is shown in Fig. 1. The characteristic superradiant pulse^"*'^^ 
is presented in Fig. 2. 

A series of numerical simulations were accomplished^^ which could be consid- 
ered as computer experiments modelling the behaviour of spins in conditions typical 
of experimental observations. These computer simulations confirmed the existence 
of both pure as well as triggered spin superradiance, being in good agreement with 
experiments. 
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3 Basic Model 



Since the phenomenon of spin superradiance exists in nature, it is necessary to 
develop its theoretical description. The simplest idea would be to invoke for this 
purpose the Bloch equations. However, these are designed so that they treat the 
magnetization of a sample as a macroscopic vector, which presupposes the existing 
coherence from the very beginning. If at the initial time coherence of spins is absent, 
the Bloch equations would never display it. Therefore, these equations can describe 
only triggered spin superradiance, when coherence is thrust upon spins at the initial 
time by assuming the presence of nonzero transverse magnetization. But pure spin 
superradiance, being a self-organized coherent process cannot in principle be treated 
by the Bloch equations. A discussion of this problem as well as the related references 
can be found in a review^^. In order to be able to describe all possible regimes of spin 
dynamics, it is necessary to resort to microscopic models. A theory of nuclear spin 
superradiance, based on realistic Hamiltonians typical of nuclear spin assemblies in 
condensed matter ' ' has been previously developed 

A system of nuclear spins, enumerated by an index i = 1, 2, . . . , A^, is charac- 
terized by their spin operators Ij. The Hamiltonian of nuclear spins in a solid can 
be written as a sum ^ 

i i^j 

of the Zeeman term and of a part due to many-body spin interactions. The Zeeman 
term contains 

H, = -h^r^B ■ ii , (2) 

where B is the total magnetic field and 7„ is the nuclear magnetogyric ratio. The 
total magnetic field 

B = Boe, + {Bi + H)e,, (3) 

consists of a static magnetic field Bq along the z-axis and of a transverse field being 
a sum of an external field B^ and of a field H produced by an electric coil which the 
sample is immersed into. The corresponding orientation of the coordinate axes with 
respect to the sample is explained in Fig. 3. The static magnetic field is directed so 
that 

InBo < . (4) 

In general, the nuclear magnetogyric ratio can be positive as well as negative. If 
it is positive, then inequality (4) tells that the static field is negative. In this case, 
the initial polarization of the sample nuclei is called inverted if it is positive, that is 
directed against the static magnetic field. 

The basic force acting between nuclear spins is due to the dipolar interactions 
which may be described by the Hamiltonian 

Hij-j:Dfiri!, (5) 

a/3 
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where the upper indices a and P denote the coordinate components a,P — x, y, z, 
while 



D?f - ^ {^c, - 3<,<.) (6) 
is the dipolar tensor, in which 



II ^ij 

fij 



The dipolar tensor satisfies the equalities 



of which the first is exact and the second is asymptotically valid for a macroscopic 
sample. 

For practical purpose, it is convenient to work with the ladder spin operators 
if" = I J ± /J, called the raising and lowering spin operators, respectively. Then, 
with the help of the notation 

c^J ^ I {D!f - loy;) , (8) 

the Zeeman term (2), for the total field (3), becomes 

= -n^nBoit - \ h^n {B, + H) (/+ + Jr) (9) 
and the dipolar interaction (5) takes the form 

+ 2c,fii^ + 2dl^iri^ . (10) 

Equations (9) and (10) will be used below to derive the evolution equations. 

The electric circuit, coupled with the spin sample, is characterized by resistance 
R, inductance L, and capacity C. The coil, surrounding the sample, has n turns of 
cross-section area Ac over a length I. The electric current j of the circuit is given 
by the Kirchhoff equation 

^1 + ^^ + ^ fmdt'-E- §, (11) 
at o JO at 

in which E is an electromotive force, if any, and $ is a magnetic flux 

$ = — nAcvM^ (12) 
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formed by the x-component of the magnetization density 

^^ = ^E<^r>- (13) 

i 

Here the brackets < . . . > imply statistical averaging. The fiUing factor 77 is approx- 
imately equal to 77 V/Vc, where V is the sample volume, while Vc = AJ is the coil 
volume. 

The electric current in the circuit is formed by moving transverse spins. In its 
turn, the current, circulating over the coil, produces a feedback magnetic field 

H = 47r^ . (14) 
cl 

Hence, the Kirchhoff equation (11) can be rewritten for the feedback field (14). To 
this end, it is useful to employ the notation for the natural circuit frequency 



and for the circuit ringing time 

~7 V~2L) 



rA f7-4l- (16) 



The related circuit damping 



7=^ fQ-^l (17) 



2Q V R 



is connected with the quality factor Q. Employing these notations, together with 
that for the reduced electromotive force 

cE 

h^^. (18) 
one obtains from the Kirchhoff equation (11) the result 

^ + 2'yH + uo^ f Hit') dt' = 7/1 - 47rr/ ^ (19) 
dt Jo dt 

for the feedback magnetic field (14). 

The feedback equation (19) can be presented in a form that extremely useful to 
exploit for analysing the evolution equations . By envolving the method of Laplace 
transforms and introducing the transfer function 

G{t) = (^cosut - ^ sinut^ e"'^* , 

with UJ = v^cu^ — 7^, one may present eq. (19) as the integral 

Git - t') hhit') - AT:r]MJt')\ dt' , (20) 

Jo L J 

where the dot over means, as usual, time derivative. 
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4 Evolution Equations 



The equations of motion for the nuclear spin operators are given by the correspond- 
ing Heisenberg equations. Prom these one aims at deriving the equations for the 
following statistical averages. The function 

u=j <ir> (21) 

defines the rotation of transverse spin components. The degree of coherence in this 
rotation is described by 

w = ^ <i+><ir>^ \u\'' . (22) 

And the average 

s^-j <!!> (23) 

gives the longitudinal spin polarization, which is analogous to the population differ- 
ence in optic systems. Keeping in mind that wavelengths at radio-frequencies are 
much larger than mean distances between spins and usually even essentially larger 
than linear sizes of a sample, one may employ the uniform approximation, assuming 
that the functions (21) to (23) do not depend on spatial variables. 

In order to obtain the evolution equations for the quantities (21) to (23), one 
averages the Heisenberg equations for the corresponding spin operators. In doing 
so, one encounters the known problem of the resulting equations being not closed 
but containing, in addition to u, w, and s, also binary averages of spin operators. 
The standard way of closing these equations would be by employing the mean-field 
decoupling, when the binary average < >, with i ^ j, is factorized onto the 
product < ff >< ij >. Such a decoupling, however, ends with the Bloch-type 
equations having the same deficiency of never exhibiting coherence if at the initial 
time w{0) = 0. This is because the mean-field decoupling completely neglects spin 
correlations leading to local spin fluctuations. Taking account of such fluctuations 
is a necessary prerequisit for describing pure spin superradiance. 

Another possibility of closing the hierarchy of spin evolution equations could 
be by writing additional equations for the binary correlators < /"/f > and by 
decoupling only the higher-order correlators, retaining untouched the binary ones. 
Unfortunately, this drastically increases the number of equations and so complicates 
the problem that it becomes practically untreatable even for equilibrium magnets^^. 
The more so for nonequilibrium systems presented by nonlinear differential equa- 
tions. 

To render the set of spin evolution equations closed so that to keep the problem 

treatable and, at the same time, to take into account local spin fluctuations, one can 

19—21 24 

apply the method of stochastic decoupling ' . The main idea of the latter is to 

separate in the cvohition equations the terms describing long-range correlations from 
the terms related to short-range fluctuations. Employing the method of restricted 
traces, one may define two types of statistical averages, one incorporating only 
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long-range correlations and another one involving short-range fields that are treated 
as stochastic variables. In what follows, the first type of averaging is denoted by 
the single angle brackets < ... >, while the second type, by the double brackets 
•C . . . Then, the binary spin correlators are presented as 

<irif>^<ir><if> (i^j): (24) 

where the averaging does not touch stochastic variables which enter the evolution 
equations in the following linear combinations: 

i^-lHil - 2?^^,/; - 2q,/;) , (25) 

with the coefficients defined in eq. (8). Stochastic variables (25) describe local 
fluctuating fields which act on neighbour spins forcing them to move. The exis- 
tence of such fluctuating flelds should lead to the appearance of dipolar dynamic 
broadening . To completely deflne the problem, it is necessary to prescribe the 
rules of calculating stochastic averages over the random variables (25). This can be 
done by considering the stochastic variables (25) as describing a white-noise process 
characterized by the stochastic averages 

< Ut) > =< > = , < UtW) > = 2r3(5(t - 1') , 

« umt') » = « mi{t') » = , « c{t)m » = ^v^sit - o , (26) 

in which is the width of dynamic broadening due to local dipole interactions 
of nuclei. This is a sort of inhomogeneous broadening existing additionally to the 
homogeneous broadening yielding the longitudinal, Fi = Tj"^, and transverse, r2 = 
T2"^, relaxation widths. 

To obtain the resulting evolution equations in a compact form, it is convenient 
to introduce the notation 

/=-i7„(Si + i/)+e (27) 

for an effective force acting on a spin. The transverse external magnetic fleld Bi 
may, in general, contain a static as well as an alternating term, as 

Bi ^ hi + hi cos ujt , (28) 

where, for simplicity, the frequency of the alternating fleld is assumed to be in 
resonance with the frequency of the electric circuit. Recall that the magnetic fleld 
produced by the coil, is caused by a feedback fleld and, possibly, by a fleld of an 
electromotive force, which may be presented as 

hit) = ho cos ujt . (29) 
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Finally, the NMR frequency is denoted by 



UjQ = I'fnBol ■ (30) 
Then the evolution equations for the functions (21) to (23) acquire the form 

— = -i{uJo + 6 - i^2)u + fs , (31) 
^^-2r2W+{u*f + ru)s, (32) 

^ = -^(«7 + r«)-ri(.-C), (33) 

where ( is an average stationary value for a z-component of a spin. Generally, 
—1<C<1. When there is no external pumping, then ( = —1. In the presence of 
pumping, e.g. by means of dynamic nuclear polarization, the pumping parameter 
C > — 1 and can reach the value C = 1- Equations (31) to (33) compose a nonhn- 
ear system of stochastic differential equations describing all dynamic properties of 
nuclear spins. 

The system of cvohition equations (31) to (33) looks yet rather complicated. For- 
tunately, there are several small parameters that allow one to simplify the consider- 
ation by employing the scale separation approach}'^~'^^''^'^ which is a generalization 
of the Krylov-Bogolubov averaging technique to stochastic and partial differential 
equations. The existing small parameters are connected with the relatively small 
values of the relaxation widths Fi and F2 as compared to the NMR frequency (30), 
so that 

— < 1 , — < 1 . (34) 

Uo UJQ 

Clearly, the quantity 

Fo = nr]pn-flhl (^pn = — ^ (35) 
is also small with respect to o^o, as well as the dynamic broadening width F3, that is 

« 1 , « 1 . (36) 

All amplitudes of the applied transverse fields are assumed to the small too, which 
implies that the values 

l/Q = 7„/l0 , J^l = Inhl , 1^2 = 7n^2 (37) 

satisfy the inequalities 

^ « 1 , ^ « 1 , M « 1 . (38) 
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The electric circuit, coupled to the spin system, is supposed to be of good quality, 
having a high quality factor, which means that the ringing width is small compared 
to the natural circuit frequency, 

^ < 1 (Q > 1) . (39) 

UJ 

The last, though not the least, is that the electric circuit has to be tuned to the 
NMR frequency (30), hence the resonance condition for the detuning A must be 
valid: 

^ < 1 (A = a; - u;o) ■ (40) 

In this way, the circuit plays the role of a resonator. 

The existence of the small parameters makes it possible, by analysing the right- 
hand sides of the evolution equations (31) to (33), to conclude that the function u 
has to be classified as fast, as compared to the slow functions w and s. Conversely, 
the functions w and s arc quasi-invariants with respect to u. In addition, one may 
derive an explicit expression for the resonator magnetic field H by iterating its 
integral representation with the solution of eq. (31), where only the main term in 
the right-hand side is retained. This iteration gives 

7„if — ia{u — u*) -\- (3 cos ut , (41) 

where the first term, with the coupling function 



a = 



7 



(l - e-^*) , (42) 



is caused by the feedback coupling of spins with the resonator, and the second term, 
with the coupling 

/5 ^ ^ (l - e-*) , (43) 

is due to an electromotive force. The time dependence of the coupling functions 
(42) and (43) describes the retardation in the resonator action on spins. Expressions 
(42) and (43) are written here for the case of small detuning, such that |A| <C 7. 
The latter assumption is not principal but is employed solely for the simplification 
of formulas. Relation (41) holds true for any detuning satisfying the resonance 
condition (40). 

Substituting relation (41) in eqs. (31) to (33) and using the notation 

/i = -ii/i -i(i/2 + /3)cosu;t-F^ , (44) 

one comes to the evolution equations 
diu 

— - = - [i{ujQ + ^0) + r2 - as] M + fis - asu* , (45) 
at 

^ = -2(12 - c^s)w + {u*h + riu)s - as [{u*f + u"] , (46) 
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^ = - ^ {u*fi + flu) -T,{s-0 + \a [{u*f + u'] . (47) 

Following further the scale separation approach, eq. (45) for the fast function u can 
be solved, with w and s being quasi-invariants. The solution obtained is substituted 
in the right-hand sides of eqs. (46) and (47), which are averaged over the period 
2-11 /uq of fast oscilaltions as well as over stochastic variables defined in eq. (25). 
Then, introducing the effective attenuation width 

^'^^ = ^'~''JI' M — + 4r l^"' ^^^^ 

where 

r = Ts - as + Tg , 
reduces eqs. (46) and (47) for slow functions to the evolution equations 

^^-2{T2-as)w + 2T,ffs\ (49) 
ds 

_ = -aw- r^ffs - r,is - c) . (50) 

In the form of the widths (48), it is again assumed that the detuning is small, 
such that |A| -C I7I -C ujq, which is not principal but just simplifies combersome 
expressions. 

The evolution equations (49) and (50) are the main equations describing the 
dynamics of strongly nonequilibrium nuclear spins. The analysis of these equations 
makes it possible to study various regimes of spin motion. 



5 Nyquist Noise 

The problem of principal importance is: What is the origin of pure spin superra- 
diance? In other words, what initiates the motion of spins when no coherence is 
thrust upon the system ai t — and there are no external fields pushing spins in 
the transverse direction? 

Keeping in mind the analogy with atomic assemblies, one may remember that 
in an inverted system of atoms the relaxation process begins with atomic sponta- 
neous radiation, which is a quantum process. After the seed radiation field appears 
in the system, atomic correlations start arizing through the interatomic photon 
exchange. As soon as these correlations become sufficiently intensive, coherence de- 
velops. Then, the quantum stage of spontaneous emission changes for the coherent 
stage, when atoms are correlated and emit coherently, which results in superradi- 
ance. 

The collectivization of spins can be produced by means of the resonator feed- 
back field. But for this field to arise, the spins, first, have to start their motion. 
Spontaneous emission for spins is absent. Then, what could initiate the motion of 
spins? What mechanism would be the origin of pure spin superradiance? 
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There existed a widespread delusion that the thermal Nyquist noise of the res- 
onant electric circuit could be the initiating cause of spin rotation. At this point, 
it is necessary to stress that the role of the thermal Nyquist noise in producing a 
fluctuating torque on the magnetization was in detail discussed by Bloembergen and 
Pound in their classic paper^. They showed that the average thermal damping is 
inversely proportional to the sample volume, because of which this damping could 
be noticeable only for a single spin, but for a macroscopic sample of many spins 

^ 1 the coil thermal damping should be negligibly small. This analysis seems 
to have been completely forgotten by later workers who have proclamed the leading 
role of the coil thermal noise in initiating spin rotation. 

To explicitly illustrate the role of the thermal Nyquist noise, one has to consider 
the effective attenuation width (48), including the influence of the electromotive 
force related to the resonance mode of the thermal noisc^^^'^^. When there are no 
external transverse magnetic fields, that is when i/i = 1/2 = 0, the effective width 
(48) reads 

r.„ = r3 + ^(i-e-")(i-e-f . (51) 

The electromotive force, corresponding to the resonator thermal noise, can be writ- 
ten as 

E = Eocosu;t. (52) 
The related magnetic field, defined in eqs. (18) and (29), has the amplitude 

fc„ = . (63) 



Prom here, it follows that 



-fRN 

With the definition (37), one gets the quantity 



hi = 81T . (54) 



2 _ SPpEg 

- h^IRN ' ^^^^ 

which characterizes in eq. (51) the average thermal damping due to the Nyquist 
noise. 

Pirst of all, eq. (51) shows that the thermal damping at i = is exactly zero 
because of the temporal factors. The latter become essential only for t > T2,t. 
Therefore, the influence of the thermal noise at short times should be suppressed by 
these temporal factors. 

Moreover, the quantity (55) is inversely proportional to the number of spins. 
That is, the thermal damping is inversely proportional to A", and thus it should 
be neghgibly small for a macroscopic sample with N ^ 1. This is in a complete 
agreement with the conclusion of Bloembergen and Pound . 

To be more precise, it is possible to explicitly calculate the quantity (55) sub- 
stituting there the square amplitude of the electromotive force due to the thermal 
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noise , which is 



At radio-frequencies, the inequahty 



is vahd. As a result, eq. (56) simphfies to 



El^- -fRjjT . (58) 

TT 



Thence, eq. (55) becomes 



2 sroco-T , . 

"° = w 

Expression (51) shows that the thermal damping has to be compared to the 
dynamic broadening width Fa, which can be of order of When the inequality 

% « 1 (60) 

i 2 

holds true, the thermal damping can be safely neglected, playing no role in spin 
relaxation, as compared to the dynamic broadening. 

To estimate the factor (59) and, respectively, to check the inequality (60), one can 

11 1 o 

accept the values typical of experiments on proton spin superradiance . Then, 
for the filling factor 77 — 0.6, proton density p„ = 4 x 10^^ cm^'^, and the proton 
magnetogyric ratio 7„ = 2.675 x 10^ s^^T^^ the linewidth (35) is Tq = 2.846 x lO'' 
Hz. At temperature T = 0.1 K, the thermal frequency is ujt = 1.309 x 10^^ Hz. 
Under the static magnetic field i?o = 1 T, the proton NMR frequency is uiq — 
2.675 X 10^ Hz. The ratio uj/ujt ~ 10~^ is small. With the transverse relaxation 
width V2 = 1.176 x 10^ Hz, one obtains 

Tl N 

For a sample volume of about 1 cm^ and N ~ 10^^, one has ~ 10"^''. This 

ratio is so much tiny that, clearly, the thermal Nyquist noise is not able to play 
any role in spin relaxation in a macroscopic sample. It is only for the number of 
spins N < 10^, when the thermal noise would be noticeable. When considering 
macroscopic samples, with ^ 10^, the resonator thermal noise is to be neglected. 



6 Incoherent Stage 

At the initial stage, the main role in starting spin relaxation is played by the dynamic 
broadening due to local spin fiuctuations. Omitting in the effective width (48) the 
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terms of second order with respect to small parameters and neglecting the resonator 
thermal noise by setting /? ^ 0, one has 

Te// = + I (l - . (61) 

This equation shows that at the very beginning of the process, when i — > 0, then 
Teff T3. That is, the dynamic broadening width Fa makes the major contribution 
to the effective relaxation width, eq. (61). 

For asymptotically small times t ^ 0, the coupling function (42) is close to zero, 
q; — > 0. Then the evolution equations (49) and (50) can be reduced to the form 

^ = -2F,^ + 2F3.% ft=-ns-n, (62) 



in which 



Fl = Fi + F3, C^^C- (63) 

i 1 



The solutions to eqs. (62) are 

w = woe-^'^'' + 2F3 /* s2(i')e-2r2(*-*') dt' , s = C + (so - C*)e"^^* , (64) 
^0 

which demonstrates the character of spin motion at short times t — > 0. At this stage, 
the motion of spins is yet completely incoherent, since, because of retardation, the 
coupling with the resonator has not yet been switched on. Assuming that Fi <^ F3, 
hence (* ^ 1, one may simplify the solution (64) as 



W ~ WqC ^^^^ 



The initial incoherent stage of spin relaxation can also be called the quantum stage, 
as far as the relaxation is caused by quantum effects of spin-spin interactions and 
there are yet no collective effects that could lead to the development of coherence. 

The duration of the incoherent quantum stage lasts till the quantum time tq, 
when the coupling function (42) reaches the value F2, so that 

a{tq) = F2 , (65) 

and when taking account of collective effects, due to the coupling with the resonator, 
becomes important. It is also easy to notice that the difference F2 — as in eq. (49) 
may change its sign at the quantum time tq, which means that the generation of 
coherent radiation would begin. Combining eqs. (42) and (65) gives for the quantum 
time 

tq^Tln^, (66) 
where the notation of the effective coupling parameter 
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is introduced. The quantum time (66) is positive only for g > 1. li g = 1, then 
tq ^ oo. This indicates that the quantum stage is finite and may change to a 
collective stage only if the coupling parameter (67) is g > 1. If the spin-resonator 
coupling is weak, and g there exists solely the quantum stage, and the collective 
stage never comes into being. For sufficiently strong coupling, eq. (66) gives 

t,^-^ (g » 1) . (68) 

The values of solutions (64) at the end of the quantum stage, i.e. at tg, can be 
estimated assuming that tq <^ T2 and <^ T3 = T^^. Then eqs. (64) yield 

W{tq) ~ Wo + 2633'^ , S{tq) ~ So , (69) 

where 

63^"^ = ^. (70) 

This is a small parameter, since Tq ~ while F2 <^ co. But, no matter how small 
^3 is, it can be principally important if Wq = 0. 



7 Transient Super radiance 

After the incoherent quantum stage, the coupling function (42) increases, switching 
on collective effects resulting in the appearance of spin superradiance^^"^"*^. An ac- 
curate description of this process can be made by solving the evolution equations 
(49) and (50) numerically. It is also possible to present an explicit analytical de- 
scription of how superradiance develops for the transient stage, when the time is 
larger than the quantum time tq but yet much smaller than the spin-lattice relax- 
ation time Ti. Then, the term with Fi in eq. (50) can be omitted. If there are no 
strong external transverse fields, then one should set i^i — > and 1/2 — > 0. Assuming 
that the coupling function (42) has reached its maximal value, one gets a ^ (7F2. 
The dynamic broadening width F3 is not larger than F2. Hence, for sufficiently large 
couphng parameter g, one may neglect F3 as compared to gT2 S> F3. Under these 
conditions, eqs. (49) and (50) read 

dw ^ ds ^ , 

- = -2r,(i-9.)», ^ = -sr.™. (71) 

Initial conditions for eqs. (71) should be taken at t = tq. 
Equations (71) can be solved exactly^^"^^ yielding 

Here Tp is a pulse width, Tp = F~^ is a pulse time, and to is a delay time. These 
parameters are the integration constants that are to be defined from initial conditions 
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w(tg) and s(tq) taken at the boundary of the transient stage, that is at the quantum 
time tg. Employing eqs. (69) results in 

= + {gT^f {wo + 2d,sl) , r, = 12(1 - gso) , TpTp = 1 , (73) 
while the delay time is 

(74) 

Solutions (72) demonstrate that the coherence intensity w{t) has the shape of a 
burst of width Tp and peaked at the delay time to- 

If the spin-resonator coupling g is sufficiently weak or the initial spin polarization 
So is sufficiently small, so that gso < 1, than the value of the delay time (74) shifts 
to the quantum region, becoming to < tq. This means that the radiation intensity 
would not have the shape of a coherent burst, but is rather a decreasing function of 
time, typical of nuclear induction. In the case of an external pulse thrusting onto 
the spins an essential initial coherence, such that g^Wo > 1, the radiation time Tp 
is small, Tp < T2, which is characteristic of collective induction. When gso > 1 and 
g'^Wo > 1, the signal of collective induction is peaked at a finite delay time to > 0. As 
is explained in the Introduction, collective induction, with coherence being induced 
by an initial external source, is in principal different from superradiance that is a 
process of spontaneously arising coherence. 

For solutions (72) to describe a superradiant burst, the requirements discussed 
in the Introduction must be satisfied. In order that the initial coherence, given by 
Wo, would not essentially influence the pulse width Tp, it should be that 

g^wo < 1 . (75) 

This guaranties a basically spontaneous character of the process. Then, the pulse 
time Tp is to be sufficiently short and the delay time finite, 

Tp < ^2 , tq<to<oo. (76) 

The second of these inequahties yields 

^so > 1 , wo + 263SI > . (77) 

Combining all conditions (75) to (77), one sees that there exist two types of tran- 
sient spin superradiance, whose classification, taking account of the smallness of the 
parameter ^3 <C 1, can be accomplished as follows: 

(i) Triggered superradiance, when 

gso>l + ^Jl- g^wo , wo 7^ . 

Here, a weak external pulse imposes an initial coherence on spins, but, being rather 
weak, this pulse plays just the role of a trigger, while the following development of 
coherence is mainly governed by internal properties. 



to^tq+^\n 



r — r 

J- p J- g 
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(ii) Pure superradiance, when 

gso > 2 



wo = . 



This is a purely self-organized process, with no coherence prescribed from external 
sources. 

In the case of pure spin superradiance, the characteristic quantities (73) and 
(74), having regard to ^3 ^ 1, become 



~ {gso - i)r2 



1 + 



.9^0 - 1 



5, 



ip — 



9S0 - 1 



1 



.9^0 - 1. 



In 




(78) 



9 2(grso - 1) 

Although the parameter ^3 is small, it cannot be neglected, since its value is crucial 
for defining the delay time to- As is seen, if ^3 0, then to ^ 00. The parameter 
^3, according to the notation (70), is proportional to the dynamic broadening width 
r3 due to local spin fluctuations. As far as the latter fluctuations are the main cause 
of starting spin relaxation, it is not surprising that the related broadening defines 
the delay time. 

From eqs. (72) it follows that the maximum of the coherent burst occurs at 
t = to, when 

wito) = (^so - , sito) = ^ . (79) 

For longer times t 3> to, the superradiant signal exponentially diminishes, and the 
spin polarization returns to an almost inverted value. 



w ~ Aw{to)e 



-So + 
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(80) 



The larger is the spin- resonator coupling g, the more complete is the spin inver- 
sion. The effect of superradiant spin inversion can be employed for the ultrafast 
repolarization of nuclear targets used in scattering experiments""^^'^^. 

The shape of a superradiant pulse, described by w{t) in eqs. (72), is in very 

11 10 

good agreement with experiments . Thus, the measured signal, presented in 

XX X2 

Fig. 2, ideally fits the calculated function w{t), as is discussed in the papers ' . 
A concrete relation between the measured intensity of the signal and w{t) will be 
explained in Sec. 10. 

It is worth emphasizing that the radiation intensity and radiation time of the 
studied pulse have the properties typical of a superradiant burst. To illustrate this, it 
is sufficient to recall that the initial inverted spin polarization so = N^/N is the ratio 
of the number of inverted spins to their total number in the sample. Inverted spins 
play the role of radiators, whose number is N+. For large spin- resonator coupling 
g ^ 1, one has w{to) ~ Sq and Tp ~ 1/so, which is clear from the consideration 
above. Therefore 



w{to) - N 



+ ' 
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which is characteristic of superradiance. 

When there is no constant pumping by means of dynamic nuclear polarization, 
there appears one dominant superradiant burst given by eqs. (72). In an inhomo- 
geneous sample, after the first burst, secondary maser oscillations can arise, whose 
accurate description requires numerical investigation ' ,ou o 

It is interesting to mention that the signals of nuclear spin echo can also be 
essentially amplified by the presence of a resonant electric circuit coupled with a 

oo 

spin system, which has been observed in experiment . 



8 Pulsing Superradiance 

Another regime of spin superradiance can be achieved if the inversion of nuclear 
spins is constantly supported by a pumping mechanism, for instance, by dynamic 
nuclear polarization. Then, a series of superradiant bursts can be produced, as is 

O 1 A 

experimentally observed . Such a regime, with a series of superradiant pulses 
that may repeatedly arise during a rather long time, can be termed pulsing superra- 
diance. The theoretical description of this regime is based on the evolution equations 
(49) and (50), which require numerical computations ' ' . However, for the late 
stage, when t ^ T2, r, and the system is close to a stationary state, it is possible to 
solve the problem analytically. 

At the late stage, when t 3> r, the coupling function (42) becomes a = gT2- As 
is explained in Section 5, the resonator thermal noise can always be neglected, that 
is /9 — > 0. When there are no strong external transverse fields, then i/i — > and 
1^2 — ^ 0. For the simplicity, the dynamic broadening width can also be neglected 
as compared to a large gr2. Then the evolution equations (49) and (50) can be 
written as a two-dimensional dynamical system 

dw ds 

with the velocity fields 

vi ^ -2r2{l - gs)w , V2^-gr2W-ri{s-C) . (82) 

In the presence of dynamic nuclear polarization, F* is a pumping rate and C > is 
a pumping parameter. 

An accurate theoretical analysis of eqs. (81) is given below. The velocity fields 
are taken in the simplified form (82). This is done for pedagogical reasons, in 
order not to become entangled into cumbersome formulas, but for emphasizing the 
principal ideas and techniques, on which such an analysis is based. Using the same 
methods, one could accomplish an analysis of the evolution equations (49) and (50), 
including the dynamic broadening and external transverse fields, which would result 
in much more intricate equations. 

An important quantity characterizing local stability of motion is the local ex- 
pansion rate"^^ 

A(t) = - Re /* TrX(t') dt' , (83) 
t Jo 
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whose definition involves the Jacobian expansion matrix X that in the considered 
case is 

dvi dvi 



OV2 OV2 

dw ds 



If, at the moment t, the rate (83) is positive, h.{t) > 0, this means that the phase 
volume of the dynamical system expands, while if A(t) < 0, then the phase volume 
contracts. The eigenvalues of the expansion matrix X are the local characteristic 
exponents 

= - ^ {rt + 212(1 - gs) ± ^J[V\ - 2r2(l - gs)f - 9,g^Vlw^ , 

whose real parts define the local Lyapunov exponents 

X^{t) = ReX^{t) . 

The signs of the latter show whether the motion at the moment t is stable or not. 

In the theory of dynamical systems , one usually considers asymptotic stabihty 
related to the limit i — > oo. The stationary, or fixed, points are given by the zeros 
of the velocity fields. The equations vi ~ V2 — 0, with the velocities (82), yield two 
stationary solutions: one fixed point is 

wl^Q, sl^C (84) 

and another is 

^; = 4^(K-i), 4 = - (85) 
g ^2 g 

As is clear, among several solutions, only those could have sense whose values are in 
the region of validity of the considered functions. This region, as follows from eqs. 
(22) and (23), is limited by the inequalities 

0<w<l, -l<s<l. 

Thence, it is evident that the fixed point (85) may have sense only if g'C ^ 1- 
When gC, — 1, both fixed points (84) and (85) coincide. On the manifold of system 

parameters, the value gC^ = 1 is termed a bifurcation point. 

Because of the existence of two stationary solutions, the limit 

A* = lim Kit) = Re lim Tr Xit) 

of the local expansion rate (83), for which A* = A"*" + A~, with the Lyapunov 
exponents 

A± = Re lim X^it) , 

possesses two different values. The local expansion rates at the first and second 
fixed points, respectively, are 

A* = _r*-2r2(l-K), k\ = -T\. (86) 
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According to the principle of minimal expansion , the smaller the local expansion 
rate, the more stable is a dynamic state. The values (86) show that when g( < 1, 
the fixed point (84) is more stable, while when > 1, the stationary solution (85) 
becomes more stable. 

The limits = limt^ao X^{t) oi the characteristic exponents, at the corre- 
sponding fixed points, are 

= xr = -2r2(i-K), 

(87) 



xi 



1±^1-8^(K-1) 



2 

The related real parts = ReX"^ define the Lyapunov exponents. Analysing eqs. 
(87) yields the following classification of the fixed points (84) and (85). 

When gd^ < 1, then the Lyapunov exponents for the stationary solution (84) are 
negative, A]*^ < 0. This tells that the fixed point (84) is a stable node. For the 
stationary solution (85), one has A^ < but A^ > 0, which means that the fixed 
point (85) is a saddle. Hence, for ^f^^ < 1, spin dynamics tends to the stationary 
solution (84). The relaxation rates are given by Xf — Xi defined in eqs. (87), from 
which it follows that the spin motion is incoherent. 

In the case g( = 1, the stationary solution is unique, since the fixed points (84) 
and (85) coincide. The related Lyapunov exponents are also the same: Aj^ = A^ < 0, 
Xi — X2 = 0. The value g( = 1 is associated with a bifurcation point. 

With the increasing pumping, when 

■p* 

one gets A^ < 0, Af > 0, while A^ < 0. This indicates that the features of the fixed 
points have been changed — the fixed point (84) is now a saddle, while that of eq. 
(85) has turned to a stable node. For realistic nuclear spins, F^ -C F2. Therefore, the 
considered region of gC is very narrow. The relaxation rates are close to A^ f» 
hence the spin motion is to be incoherent. The coherence intensity W2, given in eq. 
(85), although is not zero exactly, but, since F^ <^ F2 and gf > 1, is very small, not 
much differing from zero. 

Increasing the pumping further, so that 

^ 1 



9C>1 + 



8F2 



changes the picture qualitatively. Then the fixed point (84) continues to be a saddle, 
since A^" < 0, Aj" > 0, but the fixed point (85) turns into a stable focus, as far as its 
characteristic exponents become complex, 

X^^-^±iu;^, (88) 

with the asymptotic frequency 



^oo^^^8^(K-l)-l- (89) 
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This is the regime of pulsing spin superradiance, when there appears a long series 
of superradiant bursts lasting about the time T* = (r]|)~^. At the intermediate 
stage, this pulsing is not periodic, but becomes approximately periodic as time 
increases. The asymptotic frequency (89) defines the time interval T^o = ^ixjuj^ 
between superradiant pulses at the late stage, when the spin evolution is close to 
the stationary solution (85). This interpulse time for sufficiently strong pumping, 
with gC, ^ 1, is 

r„^^^../lSl. (90) 

The time (90), by varying the related parameters, can be varied widely. For example, 

11 1 o 

keeping in mind the parameters typical of the proton spins , one has Fq = 

2.85 X 10^ Hz, a;o = 2.675 x 10^ Hz, F2 = 1.18 x 10^ Hz. With the quahty factor 

Q = 100, the ringing width is 7 = 1.34 x 10^ Hz. The spin-resonator coupling (67) is 

g = 48. For the pumping characteristics ( = 1 and F* = 10 Hz, the interpulse time 

(90) is Too = 0.59 X 10~^ s. The number of pulses can be estimated as /Too ~ 100. 

A detailed description of the regime of pulsing spin superradiance, for arbitrary 

times t > 0, can be accomphshed with the help of numerical calculations ' ' . 

The phenomenon of pulsing spin superradiance may be employed for creating pulsing 
22 

spm masers . 

A regime, similar to pulsing spin superradiance, can also be realized without 
dynamic nuclear polarization. For this purpose, one could change in the following 
way the experimental setup used for realizing the transient spin superradiance. In 
the latter regime, as is described in Section 7, if the spin-resonator coupling g is 
sufficiently strong, a sharp superradiant burst occurs at the delay time tg. After 
the time tg + '^pj the spin polarization, according to eq. (80), becomes practically 
inverted, provided g ^ 1. If at this moment, one inverts the static magnetic field 
Bq or acts on spins by a resonant 7r-pulsc, or just turns the sample 180° about an 
axis perpendicular to Bq, then again a strongly nonequilibrium state of almost com- 
pletely inverted spins is prepared. After the time to, counted from the moment when 
the newly inverted state is created, another superradiant burst will arise. Then, one 
could again either invert the static magnetic field or invert the magnetization by 
a resonant 7r-pulse, or turn the sample to arrange a novel nonequilibrium inverted 
state of spins. After this, one more superradiant burst will appear. Such a proce- 
dure can be repeated as many times as necessary for producing a required number 
of sharp superradiant pulses. The achieved regime may be named punctuated spin 
superradiance. The principal difference of this regime from the pulsing spin superra- 
diance is the possibility of controUing the number of pulses as well as the interpulse 
time. The term punctuation here implies the possibility of varying the time intervals 
between pulses and to form groups of pulses, containing different numbers of the 
latter. In that way, it is feasible to compose a code, similar to the Morse alphabet. 
Hence, punctuated spin superradiance can be used for processing information, which 
may be employed, for instance, in quantum computers. 
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9 Hyperfine Interactions 



Real solids, in addition to nuclei, always contain electrons. It is therefore important 
to understand the influence of hyperfine spin-spin interactions between nuclei and 
electrons on nuclear spin superradiance. 

A system of nuclear and electronic spins is described by the Hamiltonian 

H — Hnuc + Hele + Hhyp , (91) 

in which the nuclear term Hnuc is given by eq. (1). The electronic spin Hamiltonian 
is 

Hele = - ^ E ' + /i7e E ^ " S. , (92) 

where Jij is an exchange interaction potential, Sj is an electron spin, and 7e is the 
electron magnetogyric ratio. The part of the Hamiltonian (91) describing hyperfine 
interactions is 

HHyp^A^S,-h + ^j:Y: <Sti^ . (93) 

Here the first term is the Fermi contact hyperfine interaction of nuclei with s- 
electrons, characterized by the energy 

^ = Y^'7e7n|V'(0)r, (94) 

with ilj{r) being the electron wave function. The second term in eq. (93) presents 
the dipolar hyperfine interactions, with 

< - 7? (^«, - 3<,<.) . (95) 

The consideration of the problem with the Hamiltonian (91) can be accomphshed 
by employing the same methods as have been detailed in the previous sections. The 
main technical difference is that the consideration becomes more cumbersome, as 
now it is necessary to deal with six evolution equations for spins, three of which are 
for nuclear spins and three other, for electronic spins ' . The seventh equation is 
the Kirchhoff equation (19), in which the magnetization density now is a sum 

M.^^Ehn <i:> -le <S! >) (96) 

of the terms due to nuclear and electronic spins. The evolution equations can be 
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again solved by invoking the scale separation approach ' ' . The most impor- 
tant conclusions resulting from the presence of hyperfine interactions are as follows. 

Local spin fiuctuations of electrons lead to the increase of the nuclear dynamic 
broadening, which results in the sum 

f 3 = Ta + r'3 (97) 
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of the widths caused by nuclear dipolar interactions, Fs, and by electron-nuclear 
hyperfine interactions, T'^. The relation between these widths is given by the ratio 

£3 ^ Pe^e 

in which pg and p„ are the electron and nuclear densities, while //g = h'jeS and 
/In = ^7n-^ a-re the electron and nuclear magnetic moments, respectively. When the 
electron and nuclear densities are close to each other, then, because fie/l^n ~ 10'^, 
the ratio (98) can be large compared to unity. Therefore, the width (97) can be 
essentially increased by the hyperfine interactions. This results in the shortening of 
the delay time to of a superradiant burst, as well as in shortening the pulse time Tp. 

In the electronic subsystem, there may appear long-range magnetic order due to 
the exchange interaction of electronic spins. If so, the NMR frequency shifts to the 
value 

ujn = uJo + — Sm^ , (99) 

where rriz is the z-projection of the average electronic magnetization normalized to 
unity, A being the parameter (94). Then, the frequency (99) will enter all formulas 
instead of uq. 

Long-range magnetic order of electrons leads to the renormalization of the cou- 
pling with the resonator according to the rule 

9^9\i + 7 — 

V Pnf^nnUJn J 

This may result because of the large ratio Pe/fJ'n ~ 10'^ in an essential increase of 
the coupling (100), as compared to g. The electronic subsystem plays the role of 
an additional resonator, which enhances the effective coupling of nuclear spins with 
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the resonant electric circuit . 



(98) 



(100) 



10 Radiation Intensity 

To complete the treatment, it is worth emphasizing once again why, actually, the 
studied phenomenon can be called spin superradiance. An ensemble of coherently 
moving nuclear spins, as is clear, generates the magnetodipole radiation with the 
total intensity 

m-^\m)\' , (101) 

where 

i 

is the total magnetization of nuclei. The quantity of interest is the radiation intensity 
averaged over fast oscillations. For the intensity (101), this yields 

m = ^^l^NMt) , (102) 
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where, as always, Hn = h'^nl is the nuclear magnetic moment. As is seen, the radia- 
tion intensity is proportional to the number of spins squared, which is characteristic 
of superradiance. 

Though coherently moving spins do produce superradiance, the radiation inten- 
sity (102) is not high. At the peak of a superradiant burst, w(to) is given by eq. 
(79). For (yf ^ 1 and Sq ~ 1, one has Wy(to) ~ 1- Then the maximal radiation 
intensity, for proton spins, with ~ 1.411 x 10~^^jr~^, — 2.675 x 10^ Hz, and 
TV ^ 10^3, is 7(to) ~ 2.5 X 10-5 W. 

Despite so weak a radiation intensity, it can be easily detected. This is because 
what is directly measured is the power of current 

P{t) - Rj\t) , (103) 

which is generated in a coil by the radiating spins. Using the relation of the induced 
current with the resonator magnetic field (14), one gets 

The field H can be found from eq. (41). Setting in this equation /3 = 0, i.e. 
neglecting thermal noise, and averaging over fast oscillations gives 



H\t) = - a\t)w{t) . 

Then, for the averaged current power (103), one finds 

Pit) = gT2inujN (l - e-^'fwit) . (104) 

This expression contains, as compared to eq. (102), an additional dependence on the 
feedback retardation. But for t ^ r, the radiation intensity (102) and the current 
power (104) differ one from another only by a numerical factor 

p(t) 3,r,A3 3gA3 

7{t) ~ IGTr^roK Stt^K ' ^ ' 

in which A = 2t:c/uj is the radiation wavelength. For an NMR frequency uj — 
2.675 X 10^ Hz, one has A = 0.71 x 10^ cm. The factor (105) can reach rather large 
values. Thus, if Q = 100 and Vc = 10 cm^, this factor is of order 10^. This is why 
even a low radiation intensity can be easily measured. 



11 Applications 

Nuclear spin superradiance is collective spontaneous radiation by nuclear spins at a 
frequency of nuclear magnetic resonance. This phenomenon is an analog of atomic 
superradiance occurring at optical frequencies. Being a novel coherent phenomenon 
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at NMR frequencies, it may find various applications, among which it is possible to 
suggest the following: 

(1) Investigation of materials characteristics 

Spin relaxation in materials is commonly employed for studying relaxation pa- 
rameters describing intrinsic properties of matter. Nuclear spin superradiance is 
accomplished by a self-organized coherent spin relaxation, which is drastically dif- 
ferent from other types of spin relaxation. Another kind of relaxation may provide 
additional information on the properties of materials. Thus, the main mechanism 
originating pure spin superradiance is the existence of local spin fluctuations caused 
by the interaction of nuclear spins with each other, by means of dipolar forces, and 
with electronic spins, through hyperfine forces. Therefore, studying specific features 
of pure spin superradiance supplies information about these local fluctuations. 

(2) Fast repolarization of targets 

The relaxation of nuclear spins, in the process of superradiance, happens much 
faster than the usual spin-dephasing time T2. If the initial spin polarization is 
sufficiently high and the spin-resonator coupling is sufficiently strong, spins, after a 
superradiant burst, change their orientation becoming almost completely inverted. 
Such an ultrafast inversion of spins can be employed for quick repolarization of 
polarized solid-state targets used is scattering experiments. 

(3) Construction of spin masers 

A superradiant spin system is a source of coherent radiation at radiofrequencies. 
Being a source of coherent radiation, it is analogous to lasers operating at optical 
frequencies. Spin masers could find applications similar to those of optical lasers. 
A spin maser can function in a transient regime, emitting a single burst typical of 
pure or triggered spin superradiance, and also it can work in a lasting regime based 
on pulsing spin superradiance. 

(4) Creation of sensitive detectors 

Spin superradiance can be triggered by a very weak external pulse of inten- 
sity corresponding to Wq. The latter depends exponentially on the delay time tg. 
Therefore, measuring the delay time of a superradiant burst provides an accurate 
evaluation of the triggering intensity. Then, triggered spin superradiance may serve 
as a sensitive mechanism for detecting weak external signals. 

(5) Methods of information processing 

In the regime of punctuated superradiance, the intervals between superradiant 
bursts as well as the number of pulses can be regulated. This makes it feasible to 
compose a kind of the Morse Code alphabet and, respectively, to develop a technique 
of processing information. Such a method of information processing could be used 
in quantum computers. 

In conclusion, it is necessary to note that the developed theory, as well as possible 
applications, are appropriate not only for nuclear spins but, in general, for any spin 
system. For instance, all this can immediately be extended to electronic spins. The 
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related phenomenon is electron spin superradiance. Other types of materials could 
be the so-called molecular magnets formed of magnetic molecules. The latter often 
have high spins, thus possessing several quantum transitions. In that case, the reso- 
nant electric circuit has to be tuned to one of the admissible transition frequencies. 
The resulting effect would be molecular spin superradiance, whose intensity could 
be much higher than that of nuclear spins. Nuclear spin superradiance is just a 
representative from a class of phenomena called spin superradiance. These phenom- 
ena, being realized with different materials, should exhibit a rich variety of specific 
features that could be not solely interesting but also useful for various applications. 
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Figure Captions 



Fig. 1. Typical experimental setup for detecting spin superradiance. The region 
surrounded by the dashed curve signifies a cryostat. Among the two coils below, 
the left one plays the role of antenna, and the right one is a part of the resonant 
electric circuit. The studied sample is shown as a dark bar inside the resonant coil. 
The upper left block is an oscilloscope, and the lower one is a plotter. 

Fig. 2. Voltage signal corresponding to a superradiant pulse as a function of 
time, measured in units of 10~^ s, for two initial spin polarizations, sq = 0.52 (lower 
curve) and sq = 0.57 (upper curve). 

Fig. 3. Orientation of the coordinate axes with respect to the sample that is 
inserted into the coil of a resonant electric circuit. 
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